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SUMMARY 


An improved transverse shear deformation theory for laminated anisotropic 
plates under bending is presented. The theory eliminates the need for an arbi- 
trarily chosen shear correction factor. For a general laminate with coupled 
bending and stretching, the constitutive equations connecting stress resultants 
with average displacements and rotations are derived. Simplified forms of 
these relations are also obtained for the special case of a symmetric laminate 
with uncoupled bending. The governing equation for this special case is 
obtained as a sixth-order equation for the normal displacement requiring pre- 
scription of the three physically natural boundary conditions along each edge. 
For the limiting case of isotropy, the present theory reduces to an improved 
version of Mindlin ' s theory . 

Numerical results are obtained from the present theory for an example of a 
laminated plate under cylindrical bending. Comparison with results from exact 
three-dimensional analysis shows that the present theory is more accurate than 
other theories of equivalent order. 


INTRODUCTION 

Classical bending theory produces errors when the ratio of the elastic 
modulus to shear modulus becomes large. Advanced cort5)osites like graphite- 
epoxy and boron-epoxy have ratios of about 25 and 45, respectively, in contrast 
to 2.6 which is typical of isotropic materials. These high ratios render clas- 
sical plate theory inaccurate for analysis of composite plates. Structures 
from advanced composites are constructed in layered form with each layer being 
orthotropic. Consequently, for analysis of composite plates, a satisf actoiry 
transverse shear deformation theory for laminated anisotropic plates is needed. 

Current shear deformation theories for laminated anisotropic plates have 
one drawback or another. The more rigorous theories (refs. 1 and 2) are cum- 
bersome because of their high order and inconvenient boundary conditions. Most 
of the simpler, sixth-order theories (refs. 3 to 5) require an arbitrary cor- 
rection to the transverse shear stiffness matrix. Cohen's sixth-order theory 
(ref. 6) does not require a correction factor but, like all other sixth-order 
theories (refs. 3 to 5), his theory is no better than classical plate theory in 
predicting stresses . 

Herein, a sixth-order bending theory for laminated anisotropic plates is 
developed that requires just the three natural boundary conditions in uncoupled 
bending problems. These three conditions are the normal moment, the twisting 
moment, and the transverse shear force. Alternatively, any one of these condi- 
tions could be replaced by prescription of the corresponding displacement 
degree of freedom, such as rotation or normal displacement. The theory does 


not require an arbitrary shear correction factor. Unlike other theories of 
equivalent order (refs. 3 to 6), the present theory gives an accurate predic- 
tion of stresses. 

The present theory uses a displacement approach like that in Mindlin's 
theoary (ref. 7) and its straightforward extensions to laminated anisotropic 
plates (refs. 3 to 5). However, in contrast to earlier theories, a special 
displacement field is used. The displacement field is chosen so that the 
transverse shear stress vanishes on the plate surfaces. This important modifi- 
cation eliminates the need for using an arbitrary shear correction factor, 
characteristic of other theories (refs. 3 to 5). 

The theory is first developed for a general unsymmetric laminate. Simpli- 
fied relations are then derived for a symmetric laminate in which bending and 
stretching are uncoupled. The results are further specialized for a symmetric 
cross-ply laminate, which is a case of classical orthotropy. For the limiting 
case of isotropy, the present theory reduces to an improved version of 
Mindlin * s theory ( ref . 7 ) . 

Cylindrical bending of a three-layered laminate is considered, as a mamer- 
ical example, to compare results from the present theory with those from the 
exact solution and other theories. 


SYMBOLS 




C 


ij 


C 


ij 


p(n) 

^ij 


"differential" bending and twisting rigidities for laminated 
anisotropic plate (i,j = 1, 2, 3) 

elements of stiffness matrix for individual lamina connecting 
stresses and strains (i=1, 2, 3, . . ., 6) 


elements of reduced stiffness matrix for individual lamina 
defining constitutive relations of plane stress type 
(i = 1, 2, 3) 


rh/2 

•'-h/2 


C . .z^ dz 
ID 


bending rigidity of isotropic plate. 


Eh^ 


12(1 - v^) 


D. . bending and twisting rigidities for laminated anisotropic plate 

(i,j = 1, 2, 3) 

E Young's modulus of isotropic plate 

Ej^,E^ elastic moduli of individual lamina in longitudinal and 
transverse directions, respectively 

G shear modulus of isotropic plate 


2 



shear moduli of individual lamina 
h laminate thickness 

k shear correction factor in Whitney and Pagano's theory 

,K 2 transverse shear stiffness coefficients for laminated anisotropic 
plate 

M ,M ,M bending and twisting moments per unit length of laminate element 
X y xy 


^x'^y'^xy membrane forces per unit length of laminate element 


Q /Q 

^x'^y 


transverse shear forces per unit length of Iciminate element 


q distributed normal load per unit area of Iciminate surface, 

positive in direction of increasing z 

s span length in numerical example 

u,v,w displacements in x-, y-, and z-directions, respectively 

U,V,W average values of displacement coitponents over thickness 

x,y,z Cartesian coordinates with z-axis oriented in thicknesswise 

direction and measured from middle plane of laminate (fig. 1) 

Px# Py average values of rotations of line normal to middle surface over 
thickness 

'^2' * * * '^11 linear differential operators 


^3^ ' ^yz ' ^xz 


I strain conponents in Cartesian coordinates 
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X 


angle of fiber orientation 


V Poisson's ratio of isotropic plate 

^Jlt ' ^tt * s ratios of individual lamina 


a ,o ,a 

X y z 


T /T ,T 
xy yz xz 


stress con^onents in Cartesian coordinates 


X transverse shear function in Reissner's theory for isotropic plates 

under bending 

Subscripts: 

Jl,t longitudinal and transverse directions of unidirectional lamina/ 

respectively 

max maximum value 

S/B stretching and bending conponentS/ respectively 

SS/BB applied to matrices/ refer to stretching and bending parts / 
respectively 

SB/BS applied to matrices/ refer to stretching-bending coupling 

Superscript : 

T transpose of matrix 


Matrix 

notation; 

c : 

matrix 

o 

column matrix 

L J 

row matrix 

CO] 

null matrix 


DEVELOPMENT OF THE THEORY 

The present theory is developed by using a displacement approach. The 
inplane displacements u and v are approximated by cubic polynomials in 
z. The out-of-plane displacement w is assumed to be constant with respect 
to z. By requiring that shear stresses vanish on the surface/ coefficients of 
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the polynomials are related. By using these relations, the displacement field 
is completely defined in terms of z, "average” displacements U, V, and W, 
and "average" rotations of a line normal to the middle surface and py. 

Then the stress distribution through the thickness is determined from the 
constitutive relations of individual layers. Next the laminate constitutive 
equations, which relate the forces and moments to displacements U, V, and 
W and rotations p^ and p^., are obtained from stresses by integration 
through the thickness. 

Finally, for the case of a symmetric laminate, the moments and transverse 
shear forces from the plate constitutive equations are substituted into the 
three plate equilibrixim equations. By eliminating all quantities except W, a 
single governing equation in terms of W is obtained. 


Unsymmetric Laminate 

The general case of an unsymmetric laminate with stretching-bending cou- 
pling is considered first. The constitutive relations for any individual 
lamina are 



( 1 ) 


(2) 


where the coordinate system is shown in figure 1. The is assumed to be 

small in comparison with other normal stresses and is neglected. Whitney and 
Pagano (ref. 4) have shown that, with this assumption, equation (1) reduces to 
the following contracted form: 



(3) 
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where 


C. . 

ij 



x4 j4 


^44 


(4) 


The strains in equation (3) are determined from an assumed distribution of dis- 
placements through the thickness* The displacements are approximated as 


2 

u = UQ(x,y) + zu^(x,y) + z U2(x,y) 

2 

V = VQ(x,y) + zv-|(x,y) + z V2(x,y) 
w = W(x,y) 


3 

+ z U 3 (x,y) 
+ z^V3(x,y) 


(5) 

(6) 
(7) 


The polynomial representation of u and v as shown in equations (5) 
and (6) is an important departure from Mindlin's theory (ref. 7) and its exten- 
sions to laminated plates (refs. 3 to 5), which assume u and v as linear 
functions of z. The assumption of a higher order polynomial in z for u 
and V in plate bending theories is not new. Lo and others (ref. 2) have 
previously used a cubic polynomial for u and v as in equations (5) 
and (6). However, in contrast to their theory, (i = 0/ 1/ 2, 3) in 

equations (5) and (6) are not independent. Instead, they are chosen so that 
the condition t^z “ ^^yz “ ^ ^ ” ±h/2 is satisfied. When these four 

conditions are satisfied, the final number of independent unknowns is five. 
Thus, the order of the theory is the same as that of earlier low-order theories 
(refs. 3 to 5). In addition, the physical conditions =0 at 

z = ±h/2 are also satisfied. It was not possible to satisfy these conditions 
in the earlier theories (refs. 3 to 5). 


Because the shear stresses x^^ and Xy^ are zero on the two laminate 
surfaces, equation (2) indicates that the shear strains also must vanish. 
Consequently, the shear strains determined from equations (5) to (7) must equal 
zero as 


9^ ^ Sw ^ ^ 3 , 2 ^ 5W 

Y = — + — = u. ± u^h + + — 

'zx 6z 9x 1 2 4 3 9x 


= 0 


( 8 ) 


9v ^ 9w ^ ^ 3 ^ 5W 

Y ± v^h t rr v^h + — 

’yz 9z dy 1 2 4 3 9y 


= 0 


(9) 
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Equations (8) and (9) lead to 


U2 = V2 = 0 


( 10 ) 


Equations (8) to (10) indicate that there are in all only five independent 
unknowns Uq/ u-| ^ Vq, , and W; however, the first four lack direct phys- 
ical interpretation and are inconvenient from the point of view of prescribing 
boundary conditions- Therefore Uq, , Vq, and v-| are replaced by four 
other quantities with physical meaning; namely, "average” inplane displacements 
U and V and "average” rotations of a line normal to the middle surface, 
which are defined as 



( 11 ) 


( 12 ) 


"Averaging" here means averaging through the thickness. Equations (11) 
and (12) are obtained from a least square approximation. Details of derivation 
of equations (11) and (12) are given in appendix A. 

By using equation (10), definitions (11) and (12) yield 


U = Uq 


(13) 

II 

> 

2 

(14) 

P = u. 
^x 1 

20 '^3 

(15) 


2 


P = V. 
^y 1 

4. 3h 

(16) 
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Solving equations (8), (9)^ (15), and (16) leads to 



From equations (3), (5), (6), and 
and the definitions of force and moment 


(17) 

(18) 

(19) 

( 20 ) 

10), the strain-displacement relations, 
stress resultants, there follows 


I T 

N = N N N 

(3^1) Lx y ^ 


I 4) 

(3x9) (9x1) 


( 21 ) 


= |m m m L 
1 ) y ^ 


M 

(3x1) 


J (|) 

(3x9) (9x1) 


( 22 ) 


where 


I = 


_SS I SB 
3x3) (3x6) 


b] 


(23) 


J = [l „ [ L ~1 

1_ BS I BbJ 

(3x3) (3x6) 


(24) 


^SS “ ^0 


(25) 


8 



SB 


- [s i '■3] 


(3x3) 


3 

(3x3) 


(26) 


^BS 


(27) 




2 

(3x3) 


I 

1 “ 4 . 

(3x3) 


.] 


(28) 



^(n) 

"'ll 

-(n) 

^12 

-(n) 

^13 

L 

n 

-(n) 

12 

-(n) 

^22 

-(n) 

23 

(3x3) 

-(n) 

*-13 

-(n) 

^23 

-(n) 

S 3 J 

5 !"’ - I 

13 J 

rh/2 

-h/2 

C . .z dz 

13 

for n 


(29a) 


(29b) 



(30) 


<•>3 " {'*’ 0 } 


(31) 


*'B 



(32) 


n 


r 

5u /9x 
n 

av /5y 


(33) 


8u /ay + av /ax 

n n 
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The matrix constitutive equations (21) and (22) are arranged so that 
matrices corresponding to stretching, bending, and stretchings bending coupling 
appear separately in the total matrix. These matrices are denoted by the s\ib- 
scripts SS, BB, and SB or BS, respectively. Also, elements of the column 
matrix <J) are arranged so that the variables signifying stretching and bending 
are grouped separately. These groups are denoted by the subscripts S and B, 
respectively. Althou^ the present formulation is for a general laminate with 
coupled stretching and bending, the arrangement of matrices allows easy extrac- 
tion of matrices for symmetric laminate relations. The symmetric laminate case 
is discussed in a subsequent section. 

By combining equations (13), (14), (17), (18), (19), and (20), the column 
matrix ^ is expressed as a function of the derivatives of U, V, W, p , 
and 6 as 

y 

4) = H (J; (34) 

(9x1) (9x11) (11x1) 


where 



(35) 


T 

6 

^B 


l^x 

ay 

5x 

av 

&y 






-JL 

dp 

■^x 


bp 


a^w 

a^w 

[dx 

dy 

&y 


bx 

ax^ 

ay^ 

ax dy 


(36) 


(37) 


«ss 

1 

0 

px4) 


px7) 

~ 0 ~ 

T 

1 

BB 

(6x4) 

1 

(6x7) 


«SS = 


1 

0 


0 0 0 
0 0 1 
1 1 0 


(39) 
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5/4 

0 

0 

-5/3h^ 

0 

0 


0 

5/4 

0 

0 

-5/3h^ 

0 


0 

0 

5/4 

0 

0 


0 

0 

5/4 

0 

0 


-5/3h^ -5/3h^ 


1/4 

0 

0 

-5/3h^ 

0 

0 


0 

1/4 

0 

0 

"5/3h^ 

0 


0 

0 

1/2 

0 

0 

-10/3h^ 


( 40 ) 


where O is a null matrix. 

Equations (21), (22), and (34) yield 


N = I H (\, 

(3x1) (3x9) (9x11) (11x1) 


(41) 


M = J H (\> 

(3x1) (3x9) (9x11) (11x1) 


(42) 


By using the strain-displacement relations and equations (5), (6), (7), 
(10), (17), (18), (19), and (20), equation (2) yields on integration through 
the thickness 


Q = 


r 





= [R] 

(2x4) 


dW/5x 




5W/5y 


(43) 


where 


R = 


[a] 

(2x4) 


5/4 

-5/3h"' 

0 

0 


0 

0 

5/4 

-5/3h^ 


1/4 

-5/3h" 

0 

0 


0 

0 

1/4 

-5/3h" 


I 

I T 
I (2x2) 
I 


(44) 
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■h/2 

*^55 


Se 

3C 

•h/2 

Se 


^66 

3C 

1 


56 


66 


dz 


(45) 


T 




^56 



(46) 


Equations (41) to (43) are the constitutive equations for a general, 
unsymmetric laminate with stretching- bending coupling. These equations define 
the force and moment stress resultants in terms of average inplane displace- 
ments U and V, normal deflection W, and average rotations of a line normal 
to the middle surface and Py. The governing equations can be derived by 

substituting from these constitutive equations into the equilibrium equations. 
Because they are lengthy and cumbersome, they are not shown herein. Derivation 
of the governing equation is carried out here only for the symmetric laminate 
case, in which there is no coupling between stretching and bending. 


Symmetric Laminate Under Bending 

For a symmetric laminate, the matrices Lg^ and L^g which couple 
stretching and bending become null matrices. Consequently, the constitutive 
equation (42) for moments simplifies to 


M 

(3x1) 


H 




BB BB ^B 
(3x6)(6x7)(7x1) 


(47) 


where the various matrices appearing on the right-hand side of equation (47) 
are defined by equations (28), (29), (37), and (40). The constitutive equa- 
tion (43) for transverse shear forces is not simplified. Further treatment of 
symmetric laminates herein is confined to derivation of relations pertaining to 
bending. The stretching relations become identical with those from classical 
lamination theory and need no separate treatment. 

The moment and force equilibrixim equations for the laminate under bending 

are 


6M /5x + 5M /5y = Q 

X ^ X 


( 48 ) 
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9M /ax + 8M /ay = Q 
xy y y 

(49) 

+ aQy/ay = -q 

(50) 


On substitution from the constitutive equations (43 ) and (47 ) , equations 
(48) to (50) yield three governing equations for three unknowns; namely, W, 
and Py. These can be reduced to a single equation in terms of W by 
eliminating p^ and py, as follows; 

Equation (43) can be rewritten as 



fdw/axl 


[p] 

(2x4) 


dW/dy 

Q„ 






where 


r -1 * 

[P] = r)\l\ 


(51) 


(52) 


The [^] and [ti] 
partitioning the 2x4 


in equation (52) are 2x2 submatrices obtained by 
matrix [R] of equation (44) as 


r ^ ^ 1 

_(2x2) I (2x2)_ 


(53) 


A column matrix Q is now defined as 


a^w 

a^w 


5Q 

X 

5Q 

X 



ax^ 


ax ay 

5x 

Sy 

5x 

by 


From equation 

(51), 

a relationship between the column matrices 

4>b 

Q is found as 



't'B = 

(7x1) 

r 

Q 

(55) 

(7x7) 

(7x1) 




where 



^11 

0 

^2 

^13 

0 

^14 

0 


0 

^22 

P2I 

0 

^23 

0 

^24 


0 

^12 

^11 

0 

^13 

0 

^14 

r = 

^21 

0 

**22 

^23 

0 

^24 

0 


1 

0 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 

0 


The elements P^j in the foregoing equation are the elements of the [P] 
matrix obtained from equation (52). By using equation (55), equation (47) is 
reduced to 


M = 0 Q 

(3x1) (3x7) (7x1) 


where 


9 = 


L 


BB 

(3x6) 


H 


BB 

(6x7) 


r 

(7x7) 


(56) 


(57) 


By using equation (56), the moment equilibrium equations (48) and (49) 
become 


''lOx - - V 

(58) 

-''42x ■ V 

(59) 


Here, A-| to Ag are linear differential operators, the coefficients of which 
are defined in terms of elements of the matrix 0. Equations (58) and (59) 
yield 
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where 




( 62 ) 


xn 


Elimination of and Q^, from equations (50), (60), and (61) results 


^4 6 

AW ^ Y. 


6^w 


• n 1 o ‘(4-i)i (4-1) i ( 6 -i)i -^(6-1)1 

x=0,1,2,..« 5 x 5 y x=0,1,2,.«* ox oy 


= ■‘I ■" ^20 ri ’"l 1 9x ay ■" ^02 7T 

ax ■' ay 


+ E *"(4-1)1 (4-?) i 

i= 0 , 1 , 2 ,... ^ byi by 


where the coefficients a and b are derived in appendix B in terms of ele~ 
ments of the matrix Q* The convention followed for the subscripts on a 
and b in equation (63) is that the two subscripts (the first one being in 
parentheses) are not to be multiplied with each other but written adjacent to 
each other. For example, ^( 6 -i)i i = 1 is ^ 5 ^* 

Equation (63) is a sixth-order governing equation for the normal displace- 
ment. It permits the three natural physical boundary conditions to be pre- 
scribed over each boundary as in Reissner's (ref. 8 ) or Mindlin's (ref. 7) 
theory for isotropic plates. Once W is determined from equation (63) and the 
prescribed boundary conditions, all the other physical quantities can be deter- 
mined in terms of W. To this end, the transverse shear forces are found 
through equations (58) and (59), moments through the matrix equation (56), and 
rotations through equation (51), 



Symmetric Cross -Ply Laminate Under Bending 

The constitutive equations for an unsymmetric laminate and a general sym- 
metric laminate were derived in matrix form so that the coefficients appearing 
in them could be calculated routinely with matrix algebra. The special case of 
a symmetric cross ply is now considered. A symmetric cross-ply laminate is 
defined as one which is comprised of only 0® and 90® layers stacked symmetri- 
cally with respect to the middle surface. In this case, there is no coupling 
between bending and twisting, and the matrix relations derived previously for a 
general symmetric laminate reduce to simple formulas. These formulas are now 
derived. 

For convenience, the x- and y-axes (fig. 1) are chosen to be oriented 
along the principal axes of the laminate. Thus 


^13 ~ ^23 ^34 ~ ° 


(64) 


C 


56 


= 0 


(65) 


and consequently some elements of the reduced stiffness matrix also vanish; 
thus 


'13 


“ *^23 “ ° 


( 66 ) 


By using equations (65) and (66), 


equations (47) and (43) reduce to 


M 

X 


11 



12 



11 


2 2 
5 W , . 9 w 

7^ 12 7T 

8y 


(67) 


M 

Y 


12 


ap 

^^•^22 



+ A 


12 


d^w 

5x^ 


+ A 


22 



( 68 ) 


M 

xy 


/5P, S?, 


33\8y 


9x 


+ 2A 




33 ax ay 


(69) 


Q 


X 




(70) 
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where 


Q = K I 
y 2 ' 


D.. 

^ 3h^ 


^.. = 15 ( 2 ) 

1 ] 4 X] 2 X] 


5 i 

rh/2 

r ^ 

- 4 j 

^55 

-h/2 

Ll - (2z/h)J 


^h/2 

r 51 

■'2 ■? J 

L/2 """" * 

[1 - (2z/h)^ 


Here can be interpreted as anisotropic plate bending rigidities. As 

shown subsequently, they become equal to the isotropic values for the limiting 
case of isotropy* The coefficients vanish for the isotropic case and 

might therefore be referred to as anisotropic "differential" plate bending 
rigidities. 

Substituting from equations (70) and (71) into equations (67) to (69) and 
using equation (50) yields 


^ p(2) ^ p(2) ^ ^ 

"" '' ‘Kjax -K^ 


„ _ _?(2) ^ _ X(2) ^ . 

y " ■'2 ^ 2 S 2 ^ 2 ^ 

ax by 




Ki /ay 


_25(2)_^^ + 

33 ax ay 23 \k^ ay ax 
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Equations (67) to (69) and (76) to (78) are two alternate forms of con- 
stitutive equations for moments for a symmetric cross-ply laminate. Equations 
(70) and (71) are the constitutive equations for transverse shear forces. 

By using equations (50 ), (76), (77), and (78), the moment equilibrium 
equations (48) and (49) can be written as 




•'' 10 ” 


- C 


1 8 x 




■i„w 


- C 


2 dy 


(79) 


(80) 


where A 3 to A.]-] are linear differential operators defined hy 


Ag = ( ) + d^ a^( )/dx^ + d^ 5 ^( )/dy^ 


(81) 


Ag = ( ) + dg a^( )/ax^ + a^( )/ay^ 


(82) 


A)o ~ ( )/ax^ + eg a^ ( )/ax ay^ 


(83) 


A)i “ ®2 )/ 9 x^ay ■'■63 &^( )/ 3 y^ 


(84) 


kJ°12 °33^ ■ K 

2 1 


(85) 
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do = 

2 33 1 


«3 - -° 33^2 


1 22 
^ ° 33 ^ ~ 


e = 

1 11 


e = 5 ^ 2 ^ + 25 < 2 ) 
®2 ^12 ^33 


e = 5 ^^^ 
3 22 


C = (D + D ) 

1 12 33 ^ 

2 


®2 k /'^12 ’’’ ° 33 ^ 

1 


Elimination of and Q from equations ( 50 ), ( 79 ), and ( 80 ) results 

in the governing equation for w as 

AgA^g( 5 W/ 5 x) + AgA^^( 9 W/ 5 y) = ^ 8 ^ 9 ^ ^ Ag( d^q/ 5 x^ ) - C 2 Aq (5 q/ 5 y ) ( 94 ) 


By substituting for the differential operators/ equation ( 94 ) becomes 
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4 4 4 

a w . ^ aw . aw 

e. — 7 + 2e« — ;; — 7 + e. 


1 4 ■ ‘■'2 2 2 • ^3 4 ^ 6 ^4®l] ^^4 2 

dx 6x 5y dy 9x ■- -’ 


6 

a w 


9x ay 


+ |d,e + e2(d + Vs ^ 

8x ay ay 


b' 


_ £. - ^ 

= q + (^1 + ^3 - + (d 2 + d^ - C 2 )— I 

ax ay 


4 4 

M [5i(^4 “ ®2^ ■" V 3 ■ VlJ~2fr 

ax '- -* ax ay 


V^4 - °2>M 

ay 


(95) 


Equation (95) has no derivative involving odd number of differentiations 
with respect to x or y. In contrast, these derivatives appear in the 
governing equation (63) for a general symmetric laminate. 


Reduction to Isotropic Plate Relations 
For the special case of isotropy 


"'ll = S 2 = ■ '' ^ 


C 12 = VE/(1 - V ) 


11 


D 22 = = C 22 ^ = Plate bending rigidity, D 


Eh 

12(1 - v^) 


°12 = "^12 = 


''33 = ^3^ = 


(96) 


(97) 


(98) 


(99) 


(100) 
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0 


( 101 ) 


‘11 


= A 


22 


==^12 = 


K 


1 



( 102 ) 


By using equations (96) to (102), the constitutive equations (67) to (71) 
and (76) to (78) reduce to a form similar to Reissner’s (ref. 8). The only 
discrepancy involves terms of q. For the homogeneous case q = 0, the two 
theories are identical. The discrepancy for q 5 ^ 0 vanishes if the contribu- 
tion from a to strain energy is neglected in Reissner's energy formulation 
(ref. 8). As this contribution is relatively small, the discrepancy between 
the two theories can be considered to be a higher order effect and negligible. 

The discrepancy in terms involving q is a consequence of the assumption 
of constant w throu^ thickness in the present theory. Reissner’s isotropic 
plate theory is based on exact satisfaction of the equilibrium equation in the 
z-direction for all z which implies variation of w with z. 


DISCUSSION 

A sixth-order governing equation for W is obtained here for a symmetric 
laminate. This is in contrast to the Reissner (ref. 8) and Mindlin (ref. 7) 
theories for an isotropic plate which give a fourth-order equation for W 
together with an auxiliary equation of second order for a transverse shear 
function x* However, the total order is the same in both cases / thereby 
requiring prescription of the same number of boundary conditions. 

A close inspection reveals that/ for the limiting case of isotropy, the 
present theory also leads to lower order equations for W and % like the 
Reissner (ref. 8) and Mindlin (ref. 7) theories. This happens because, for 
isotropy, the differential operators Ag and Ag in equations (79) and (80) 
become identical. Thus, fewer differentiations would be required to elimi- 
nate and Q from equations (50), (79), and (80). (See derivation of 

eq. (94).) 

The and for a laminated plate are determined completely in terms 

of W as particular^ integrals of equations (79) and (80). Complementary solu- 
tions of these equations are not admissible as they violate the equilibrium 
equation (50). Recall that, in contrast/ complementary solutions are used in 
isotropic plates (refs. 7 and 8) because coitplementary solutions in the iso- 
tropic case can be chosen to satisfy the equilibrium equation (50). 


Shear Correction Factor 

The term "shear correction factor” in transverse shear deformation 
theories is usually meant to refer to an arbitrary correction applied to the 
shear stiffness previously determined. In the present theory/ there is no 
shear correction factor in this sense of the term because the transverse shear 
stiffness is explicitly determined and no correction is necessary. This 
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becomes clear by examining the limiting case of isotropy in detail* From equa- 
tions (70), (71), and (102), transverse shear forces for an isotropic plate are 
given by 



(103) 


(104) 


where G is the isotropic shear modulus. Equations (103) and (104) agree with 
the results from Reissner's theory for isotropic plates (ref. 8). 


The quantities 


3 + 


5W 

5x 


and 


aw 

9y 


represent average shear strains 


through the thickness. Thus, the factor 5/6 (0.833) in equations (103) 
and (104) can be looked upon as a correction to be applied to the transverse 
shear stiffness to account for variation of shear stress through the thickness. 
In Mindlin's theory (ref. 7), the assumption of constant shear strain led to a 
factor of unity instead of 5/6 (0.833) in equations (103) and (104). However, 
Mindlin replaced this factor of unity by an arbitrary factor, the value of 
which was adjusted so that results from the theory agreed with the exact solu- 
tion for a chosen example. He arrived at a value of *n:^/12 (0.822) for this 
factor by considering one example. By considering a second example, he 
obtained another value which depends on Poisson’s ratio v and varies from 
0.76 to 0.91 as v varies from 0 to 0.5. The first value 71^/12 (0.822) is 

close to 5/6 (0.833) obtained from Reissner’s theory. However, the manner in 
which the shear correction factor is derived in Mindlin 's approach is arbitrary 
because the value arrived at depends on the example chosen. The present theory 
is also based on a displacement formulation similar to Mindlin 's but the shear 
correction is derived in a logical way. Thus, the present theory can be looked 
upon as an improvement of Mindlin ’s theory for isotropic plates. 


Straightforward extensions of Mindlin’ s theory to laminated composite 
plates (refs. 3 to 5) have the same degree of arbitrariness as Mindlin ’s theory 
itself. For instance, Whitney and Pagano (ref. 4) arrived at different shear 
correction factors for two-layered and three- layered plates. The suggested 
procedure in their method appears to be to derive the shear correction factor 
for each set of lamination parameters by considering the known exact solution 
for a certain problem. But with the present theory, no shear correction factor 
is necessary and the transverse shear stiffness is obtained as a function of 
elastic constants and the stacking sequence without considering the exact solu- 
tion for a specific problem. The required expressions are given by equa- 
tions (74) and (75). 


22 



Interlaminar Shear Stresses 


Interlaminar shear is one of the sources of failure in laminated plates. 
Therefore, calculation of interlaminar shear stresses is important in any bend- 
ing problem. With the use of the present theory, these stresses can be calcu- 
lated as follows: A problem is considered solved if p^, and W are 

determined as functions of x and y. The displacements can be determined as 
functions of x, y, and z by the use of equations (10), (17) to (20), and 
(5) to (7). Thus, cr^, Cy, and are also determined as functions of 

X, y, and z from the constitutive relations of individual laminae. The 
interlaminar shear stresses are then determined from the following equilibrium 
equations: 



This method gives single-valued interlaminar shear stresses at the inter- 
faces. An alternate method would be to calculate and Ty^, directly from 

the displacements given by equations (5) to (7). However, this method gives 
two values for the interlaminar shear stress at each interface depending upon 
the lamina chosen. 


NUMERICAL EXAMPLE 

A numerical example is used to compare the present theory with the exist- 
ing theories. The example chosen is that of cylindrical bending of a three- 
layered, symmetric cross-ply laminate (0®/90°/0°) of high modulus graphite- 
epoxy. The layers are all taken to be of equal thickness with fibers in the 
outer layers oriented in the direction of bending (fig. 2). The layer proper- 
ties used are 
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The plate is considered to be semi-infinite with a finite span s in the 
x-direction (fig* 2) and subjected to a sinusoidal, distributed normal load of 
intensity q = ^max boundary conditions considered are 


At X = 0 and s , 

W = 8 = M =0 (105) 

•^y X 

For the case of cylindrical bending, all derivatives with respect to y 
are zero* Thus equation (95) yields 


r-t r, • 

W = W sin — 

max s 


(106) 


where 


1 + (c - d - d )(ic/s)^ + " c.)(ir/s)^ 

W = q (107) 

max ^max ~ “ 


e^(n/s)^ 1 - d^(ir/s)^ 


The transverse shear force is obtained from equation (79) as 


Q = Q cos 

X x,max s 


TUX 


(108) 


where 


x,max 


in/s) e.W - (7t/s)c.q_ 
1 max 1 nma> 

1 “ (n/s)^d. 


(109) 


Equation (70) yields the rotation 


6 as 

X 


6 =6 cos 

^x ^x,max 


1CX 

s 


( 110 ) 
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where 


3 = Q “ (ti/s)W 

f max K ^ x , max max 


( 111 ) 


From equations (80) and (71) and the boundary condition on Py 


Q = 0 P = 0 (112) 

y 

Stresses in the different laminae are now determined as follows. In pure 
bending, the terms Uq and Vq in equations (5) and (6) vanish. Conse- 
quently, the inplane displacements u and v are defined, in view of equa- 
tion (10), by 


u 


u^z + u 



(113) 


V 


v^z + v^z 


3 


(114) 


The u^ and U 2 are determined from equations (17) and (18). It follows 
from equations (19), (20), (106), (107), and (112) that v-j and v^ are zero. 
The inplane displacements are thus determined. Bending stress distribution 
through the thickness is now determined from the constitutive relations of 
different laminae. 

Exact solution for this problem based on three-dimensional elasticity 
analysis was given by Pagano (ref. 9). Figures 3 and 4 compare results from 
the present theory with those from the exact solution and the theory of Whitney 
and Pagano (ref. 4). The theories presented in references 3 to 5 are all a 
simplified set requiring an arbitrary shear correction factor. In the results 
presented in these references, the value of the shear correction factor is 
adjusted so that the results come close to the exact solution. Therefore a 
comparison of results from all these theories could be misleading. For this 
purpose, only Whitney and Pagano *s theory (ref. 4) is chosen for comparison and 
is treated as being representative of the simplified theories (refs. 3 to 5). 

Figure 3 shows a plot of the maximum deflection W^^ax center of 

the span as a function of the span-to-thickness ratio. The present theory is 

closer to the exact solution than Whitney and Pagano’ s theory (ref. 4) with a 

shear correction factor k of unity. On the other hand, Whitney and Pagano ’s 
theory shows better correlation if k is shown as 2/3. However, the factor 
2/3 was arrived at by Whitney and Pagano hy a trial-and-error procedure. There 
are two disadvantages in Whitney and Pagano ’s method. First, the trial-and- 
error procedure is to be repeated if the lamination parameters are different. 
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There is no single value of k which holds good for all lamination parameters. 
(For example, Whitney and Pagano arrived at another value for k, namely, 5/6, 
for a two-layered plate. ) Secondly, the shear correction factor arrived at by 
this procedure is problem- dependent and it is not sure whether the value for 
k so derived is valid for a problem other than cylindrical bending. Recall 
that Mindlin (ref. 7) arrived at different values for k by considering dif- 
ferent problems. 

Figure 4 presents the bending stress distribution through the thick- 

ness for a laminate with a span-to-thickness ratio s/h of 4. The exact solu- 
tion of Pagano (ref. 9) and results from Whitney and Pagano 's theory and the 
present theory are presented in this figure. Results from the classical lami- 
nated plate theory are also included for comparison. The exact solution and 
the present theory show considerable deviation from the classical theory. 
Furthermore, the exact solution and the present theory are in good agreement 
for most of the thickness. Note that in the middle layer, the 90® lamina, the 
stresses cr^ are extremely small and all theories predict near zero values. 

Figure 4 shows that Whitney and Pagano' s theory (ref. 4) yields the same 
stress distribution as the classical laminated plate theory irrespective of the 
value of the shear correction factor used. This is true of all current sixth- 
order theories for laminated anisotropic plates (refs. 3 to 6). The fact that 
these theories predict stresses no different from the classical theory is a 
severe limitation. Even Cohen's theory (ref. 6), which predicts W with good 
accuracy without a shear correction factor, has this drawback. A careful 
inspection reveals that this drawback in current sixth- order bending theories 
is a consequence of the assumption of linearity of £y/ and Yjqt with 

respect to z. The present theory allows for a more general variation of these 
strains with respect to z and does not have this drawback. 


CONCLUDING REMARKS 

A shear deformation theory for laminated anisotropic plates is devel- 
oped. In the case of uncoupled bending, the present theory is one of sixth 
order and requires just three natural boundary conditions. Most of the current 
sixth-order theories require an arbitrary shear correction factor and all of 
them have the drawback that their stress prediction is highly inaccurate. The 
present theory does not have either of these drawbacks. 

The present theory is not presented as an improvement over the current 
higher order theories. Surely, they should be more accurate but they require 
prescription of inconvenient boundary conditions. From the engineering point 
of view, it is difficult to prescribe anything other than the three natural 
boundary conditions involving moments and forces or the corresponding rotations 
and displacements. Consequently, a sixth-order bending theory which requires 
the three boundary conditions is desirable. This paper was aimed at developing 
the best possible theory restricting the order of the theory to six. 

The theory is developed for three cases. First, the formulation is 
carried out for an unsymmetric laminate. Here, the constitutive equations con- 
necting moments and forces to average displacements and rotations are derived. 
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Second, simplified forms of these constitutive equations are derived for a 
symmetric laminate. In both these cases, the various relations are derived in 
matrix form to avoid tedious algebra and also to facilitate routine computation 
on a computer. Last, further simplified relations are derived for a symmetric 
cross ply which is a case of classical orthotropy. In this case, it was found 
that matrix formulation could be dispensed with and the various relations are 
obtained in the form of simple formulas. 

For the limiting case of isotropy, the present theory reduces to an 
improved version of Mindlin's theory. The shear correction factor of 5/6 in 
this case is obtained from the theory rather than with a specific example as in 
Mindlin ' s approach . 

Of particular interest should be the method of confuting the interlaminar 
shear stress suggested in this paper. The use of equilibrium equations leads 
to single-valued shear stresses unlike straightforward confutation from 
displacements . 

The accuracy of the present theory is demonstrated by considering a 
numerical example of cylindrical bending of a three- layered plate. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
June 30, 1981 
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APPENDIX A 


DEFINITION OF "AVERAGE” VALUES OF DISPLACEMENTS AND ROTATIONS 

Let the displacements u and v be arbitrary functions of z. Let U 
and V represent average values of displacements and and py represent 

average values of rotations* "Averaging" here means averaging through the 
thickness. Essentially, the deformed shape of a line normal to the middle 
surface is sought to be approximated by a straight line so that its deviation 
from the true deformed shape is minimum. To this end, the displacements u 
and V are expressed as U + p^z and V + PyZ, respectively. Then, the devi- 
ation from the true deformed shape is represented hy the following integrals of 
squares of errors in determining u and v: 


E 

u 



U - p z) 


2 


dz 


E 

V 



(v 



dz 


The best least square approximation is that which satisfies the conditions 


5E 5E 5E 9E 

u __ u _ y _ y 

au ” 6p 5V 5P 

•^x y 


(A1 ) 


Equations (Al ) yield the definitions given by equations (11) and (12). The 
definitions given by equations (11) and (12) are identical with those obtained 
by Timoshenko and Woinowsky-Krieger (ref. 10) from considerations of work done 
by the forces and moments. However, the derivation in reference 10 is valid 
only if the stress varies linearly through the thickness. The stress variation 
for a laminated plate can, at the best, be only piecewise linear. In such a 
case, it would be necessary to resort to a mathematical definition as given 
herein. 
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APPENDIX B 


DERIVATION OF COEFFICIENTS IN THE GOVERNING EQUATION 
FOR A SYMMETRIC LAMINATE 


in* 


Elimination of Q and Q from equations (50), (60), and (61) results 
^ y 


+ A2Ag) + ^(A^Ag + A3A^^W + A^q = 0 


By substituting for the differential operators A-j to in terms of ele- 
ments of the matrix [0] , the following e3q>ressions for coefficients in the 
governing equation (63) are obtained: 


«60 - ® 16«31 - ® 11®36 


*51 ■ » 11»34 ♦ « 17«31 * » 16‘»21 ^ » 33 > ‘ » 13®36 ' ® 14»31 ‘ » 11<»26 * ® 37 > 


*42 ■ Sl®37 * ®11<®24 * ®35> * ®16*®23 * '’32’ * ^4**13 * ®3l’ 

* <®17 ® 36><®21 ® 33 > - ^ 1«27 ' ® 14<®21 « 33 > ’ ®31 <«15 ^ » 34 > 


- ® 36‘®12 ® 33 > - ‘^3 * « 31>‘®26 * « 37 > 


•33 ■ ® 11®25 * '* 16^22 * ® 34‘®12 * ® 33 ’ * ® 37<®21 ® 33 ' 


+ (6^3 + ®3 i >(024 ®35* * ^®17 ®36*^®23 ®32* " ®31 ®35 ” ®32®36 


» 14'»23 * ® 32 ’ - » 27<»13 * ® 31 > ' <»12 » 33>'®26 ® 37 > 


«15 * ® 34>'®21 » 33 > 
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® 32®34 ^ ® 22^®17 ^ ^ 36 > ^ ® 25^®13 ^ ® 31 > ^ ® 37<®23 ^ ® 32 ^ 


+ ( 0^2 ■'■ ® 33^^®24 '*’ ® 35 ^ “ ® 14®22 " ® 27^®12 '‘’ ®33 ^ 


® 32^®26 ‘’’^ 37 ^ ” ® 35^®21 ■*■ ®33 ^ " ^®15 ''' ® 34^^®23 ®32 ^ 


® 24®32 '*’ ® 25^®12 ''’ ®33 ^ ® 22^®37 ” ®15 " ®34 ^ " ® 23®35 " ® 27®32 


® 25®32 " ® 22®35 


'11 


®13 ^ 203 , 


®12 ®21 2033 


®23 2032 


"22 


®14 ^ ®36 


®15 ^ ®26 ®34 ®37 


®27 ^ ®35 


0-l/c0'3>1 ■ 0<>10-»C 
16 34 14 36 
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Figure 3.- Maximum deflection as function of thickness parameter. 


Nondimensional bending stress 
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